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道 流 条 件 下 仿生 鱼 自主 游 动 行为 研究 


周 建 军 , 袁 显 宝 , 毛 璋 亮 , 肖 仁政 
(三 峡 大 学 机 械 与 动力 学 院 ,443000 宜昌 ) 


摘 要 :深入 开展 鱼 类 自主 游 动机 理 研究 ,对 鱼 类 运动 行为 以 及 仿生 技术 研究 都 具有 重要 的 意义 。 
以 某 型 仿生 鱼 为 研究 对 象 建立 了 鱼 体 的 简化 模型 和 鱼 体 摆动 的 计算 模型 ,采用 动 网 格 技术 开发 了 
控制 鱼 体 摆 动 的 UDF 程序 ,利用 Fluent 软件 对 仿生 鱼 的 自主 游 动 行为 进行 了 仿真 模拟 ,研究 了 和 鱼 
类 自主 游 动 的 机 理 以 及 在 逆流 情况 下 不 同 来 流 方向 时 的 受 力 情况 。 通 过 计算 发 现 , 逆 流 时 水 流 方 
向 对 仿生 鱼 的 自主 游 动 有 很 大 的 影响 , 当 水 流 方 向 和 鱼 游 动 方 向 成 180°, 鱼 向 前 游 动 的 阻力 最 小 ， 
侧 向 分 量 的 存在 会 减 小 鱼 向 前 运动 的 驱动 力 。 
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Numerical analysis about hydrodynamics characteristics 
for autonomous swimming of bionic fish in counter flow 


ZHOU Jianjun, YUAN Xianbao,MAO Zhangliang,XIAO Renzheng 


( College of Mechanical & Power Engineering of China Three Corges University ,443000 Yichang ,China ) 


人 Abstract :It is very important to study the swimming rules of fish for biomimetic robotic and fish behavior 
-iesearch. This article established the simplified model and calculation model about fish oscillation based on 
a kind of bionic fish and self-propelled swimming. The calculation used dynamic grid technique and devel- 
oped UDF program for fish oscillating control ,analyzed the self-propelled swimming mechanism of fish with 
software Fluent, and also analyzed the influence of inflow angle of fluid on fish swimming behavior in 
counter flow. It can be found that inflow angle has a larger effect on fish swimming in counter flow ,the re- 
sistance of fish swimming is smallest when the angle between inflow and fish swimming direction is 180. 

Lateral component of inflow can decrease the driving force of autonomous swimming. 


Key words :self-propelled swimming;bionic movement; dynamic mesh; CFD simulation 


在 过 去 几 十 年 中 ,人 们 在 研究 自然 界 动物 的 推 ” 动 性 , 且 运 动 效 率 高 .隐蔽 性 好 等 优点 ,吸引 了 大 量 
进 机 制 方面 做 了 大 量 的 工作 ,这 对 于 解决 空气 动力 ”的 研究 者 对 其 运动 机 理 进 行 研究 。 
学 中 的 一 些 问题 以 及 在 仿生 机 器 设计 方面 都 具有 重 部 栋 伟 等 ”基于 浸入 边界 法 的 流 固 耦 合 建 模 思 
要 的 意义 。 其 中 鱼 类 由 于 在 水 中 具有 非凡 的 机 。” 想 ,以 鱼 体 肌 肉 提供 的 主动 力 为 原 动 ,建立 了 鱼 体内 
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力 \ 鱼 体 运动 和 外 界 流体 耦合 作用 的 自主 游 动 柔性 
鱼 模型 ,对 不 同 长 度 尾 鱼 , 不 同 弹性 模 量 鱼 体 的 自主 
直线 巡游 进行 了 模拟 。 严 慧 云 等 中 利用 Fluent 软 
件 , 采 用 动 网 格 技术 对 仿生 鱼 的 自主 游 动 进行 了 模 
拟 分 析 , 结 果 表明 在 自主 巡游 阶段 ,仿生 鱼 通过 摆动 
鱼 体 ,在 鱼 体 两 侧 流 场 产 生 高 压 和 低压 区 ,依靠 流体 
的 正 压 梯度 获得 前 进 的 推动 力 ,在 C 形 快速 起 动 阶 
段 , 仿 生 鱼 通过 快速 大 幅度 的 回 摆 获 得 加 速度 ,从 而 
实现 快速 起 动 。Deng 等 “利用 混入 边界 法 ,研究 了 
单 体 鱼 和 多 体 鱼 的 自主 游 动 行为 。 王 亮 等 "采用 

自 适应 浸没 边界 方法 ,数值 模拟 了 仿生 鱼 在 不 同 自 
由 度 个 数 情况 下 的 游 动 , 对 比分 析 了 自由 度 个 数 对 
仿生 鱼 受 力 、 能 耗 、 推 进 效率 以 及 尾 涡 结 构 的 影响 。 
研究 结果 表明 ,推进 效率 随 着 自由 度 个 数 的 增加 而 
略 有 降低。Cui 等 "采用 游 动 指数 对 BCF 鱼 的 游 
动 答 为 进行 研究 ,分 析 了 抽动 频率 波形 等 参数 对 游 
动 的 影响 。 夏 全 新 等 "采用 有 限 体积 法 ,并 利用 动 
枪 技 术 对 鱼 类 波状 摆动 推进 进行 了 数值 模拟 , 研 
究 巴 零 厚 度 波纹 板 的 摆动 推进 以 及 用 NACA 标准 恤 
现 柜 拟 有 厚度 鱼 体 的 波状 摆动 。 徐 晓 峰 '” 通过 运 
用 羽 工 伪 压 缩 方法 求解 三 维 不 可 压缩 非 定常 N-S 方 
稳 , 对 人 金枪鱼 和 鲁 鱼 的 自主 游 动 以 及 鲤鱼 游 动 轨迹 
的 骤 制 问题 进行 研究 ,分 析 了 影响 鱼 类 游 动 的 因素 ， 
并 过 鱼 的 游 动 控 制 提出 了 一 些 方法 。 陈 宏 等 5] 基 
也 EC 形 转向 的 仿生 学 研究 和 简化 物理 模型 ,建立 了 
鲁 暴 机 器 人 快速 转向 的 运动 学 方程 ,分 析 体 干 中 线 
运 加 步 态 和 和 鱼 体质 心 的 运动 规律 ,提出 了 有 效 提高 
鱼 形 机 器 人 快速 转向 性 能 的 方法 。 

加 虽然 针对 鱼 类 运动 人 们 开展 了 大 量 的 研究 也 取 
得 了 一 定 的 研究 成 果 , 但 是 鱼 游 动 的 内 在 机 理 是 非 
常 复杂 的 ,不 仅 要 考虑 鱼 的 主动 控制 .与 流体 相互 灶 
合作 用 等 问题 ,还 要 考虑 到 不 同 鱼 类 的 不 同 鱼 体外 
形 及 游 动 特点 。 目 前 针对 鱼 类 运动 行为 的 研究 ,大 
部 分 都 集中 在 鱼 体 游 动 行为 的 机 理 以 及 本 体 运 动 参 
数 的 影响 "5 ,对 于 游 动 过程 中 外 部 环境 的 影响 研 
究 的 较 少 。 而 鱼 类 在 实际 的 游 动 过 程 中 外 部 的 流体 
并 不 会 处 于 理想 状态 ,流体 参数 都 是 随时 发 生变 化 
的 ,在 这 种 情况 下 研究 鱼 体 的 受 力 和 运动 情况 ,对 于 
分 析 流体 对 鱼 类 运动 的 作用 机 理 和 研究 仿生 鱼 的 运 
动 控 制 都 有 重要 的 意义 。 由 于 鱼 在 水 中 的 游 动 可 以 
近似 看 作 是 在 平面 上 的 二 维 运动 ,而 鱼 游 动 过 程 中 
的 流 场 分 布 却 是 三 维 的 。 因 此 ,本 研究 以 尾鳍 游 动 
的 鱼 为 研究 对 象 建立 了 仿生 鱼 的 三 维 分 析 模型 , 研 
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究 了 在 二 维 运动 规律 下 仿生 鱼 的 自主 游 动 以 及 来 流 
方向 对 仿生 鱼 自主 游 动 行 为 的 影响 。 


1 物理 模型 


参考 常见 尾 鱼 推进 模 式 鱼 的 实体 模型 ,忽略 胸 
鳍 \ 腹 鳝 .背鳍 等 结构 ,建立 仿生 鱼 的 三 维 分 析 模 型 ， 
见 图 1。 鱼 体外 形 为 半 椭 圆 + 抛物 线形 流 线 , 鱼 体 
三 维 尺寸 为 :x xyxz=0.25 m x0.07 m x0.03 m, 图 
中 :wv 为 仿生 鱼 自主 游 动 速度 ;uw 为 流体 流动 速度 ;9 
为 仿生 鱼 游 动 方向 和 来 流 方向 的 夹 角 。 


多 1 鱼 体 模型 
Fig.1 Model of fish 
在 本 研究 的 模拟 过 程 中 只 考虑 仿生 鱼 在 二 维 运 
动 规律 下 的 自主 游 动 ,这 种 运动 规律 能 够 反映 鱼 体 
在 正常 游 动情 况 下 的 运动 行为 ,仿生 鱼 体 绕 流 的 控 
制 方程 ”为 
V =0 


1 
tu Vu Vp tv vu 局 


其 中 :uw 为 流体 速度 ;p 为 流体 压力 ;p 为 流体 密度 ;v 
为 流体 的 动力 学 黏 性 系数 。 依 据 仿 生 鱼 自主 游 动 状 
态 ,可 以 得 出 其 受 力 控制 方程 为 
Wh (2) 
其 中 :m 为 鱼 体质 量 ;V 为 鱼 体 游 动 速度 ;F 为 鱼 体 
在 速度 方向 上 受到 的 合力 ;1 为 力矩 ;L 为 动量 矩 。 
了 = DP mr xm (3) 


其 中 :m; 是 第 i 个 单元 的 质量 ;7; 是 单元 i 的 中 心 ;v; 
是 鱼 体 坐标 中 的 位 置 变 化 速率 ,v; 由 两 部 分 组 成 , 即 
DV; = + ;XT (4) 
其 中 :vw 由 摆动 规律 决定 ;w, 是 单元 的 角速度 。 
因此 
ES Dmir; x va + Dmr x ow, (5) 
为 了 使 式 (2) 有 唯一 解 ,必然 有 


Wi = = = (6) 
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因此 
du 
A 二 下 ， 
< 人 ( )+ 二 ?)=M (7) 
dt 2 Xv di OO mr 加 
当 w 和 w 确定 后 , 鱼 体 新 的 位 置 和 新 的 迎 角 就 
可 以 确定 。 
= (8) 


其 中 :9=0,i+0, + 0.k,0,、9, 和 0, 分 别 表示 在 y-z、 
x-z.x-y 平面 鱼 体 中 线 和 坐标 轴 y\z\x 正方 向 的 
夹 角 。 

为 了 减 小 计算 边界 对 仿生 鱼 自主 游 动 行为 的 影 
响 , 本 研究 取 鱼 体 长 度 为 1, 计 算 区 域 设 置 为 71 x31。 
值 林 在 高 度 上 处 于 计算 区 域 的 中 间 位 置 ,在 长 度 方 
向 下 鱼 头顶 点 位 置 距离 计算 区 域 和 人口 边界 4/, 鱼 尾 距 
离 放 算 区 域 出 口 边界 27。 计 算 区 域 设 置 如 图 2 所 示 。 


图 2 计算 区 域 

了 Fig.2 Computed region 
由 于 鱼 体 模 型 结构 的 不 规则 ,在 划分 网 格 时 为 
了 保证 网 格 质量 ,采用 四 面体 网 格 对 计算 域 进行 网 
格 蚜 分 ,同时 对 仿生 鱼 表面 进行 了 局 部 网 格 加 密 , 计 
算 虐 网 格 划分 由 Gambit 软件 完成 。 利 用 Fluent 软 
件 移动 网 格 功能 来 模拟 鱼 体 的 运动 ,仿生 鱼 运动 规 
律 的 控制 通过 编写 相应 的 UDF 程序 实现 ,参考 相关 

文献 取 仿生 鱼 轴线 运动 规律 为 
h(x,t) =a(x)sin( kx — wt) (9) 
a(xX) =cx, + Ci) ， -0. 194 <x, <0.056 (10) 
其 中 :h(%,t) 为 仿生 鱼 中 心 线 的 侧 向 位 移 ;x 为 仿生 
鱼 中 心 线 在 体 长 方向 上 的 位 置 坐标 ;a(x ) 为 波幅 包 
络 线 函数 ;= 2m/A 为 鱼 体 波 的 数量 ,文中 取 入 = 
1.4;w 为 仿生 鱼 摆 动 的 角 频 率 ,w =2mf,f 为 仿生 鱼 

摆动 的 频率 ,本 研究 取 1。 

入 口 边界 条 件 设置 为 速度 入 口 ， 速 度 大 小 为 来 
流速 度 ， 出 口 边界 设置 为 自由 出 流 。 在 分 析 仿生 鱼 
自主 游 动 时 ， 不 考虑 来 流 影响 ， 进 出 口 均 设置 为 辟 
面 边 界 条 件 ; 时 间 步 长 设置 为 0.01 s， 计 算 时 间 为 
6 s, 
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2 结果 与 讨论 


2.1 自主 游 动 行为 分 析 


仿生 鱼 在 运动 过 程 中 ,周围 流 场 分 布 会 对 仿生 
鱼 的 运动 产生 一 定 的 影响 , 流 场 分 布 对 仿生 鱼 运动 
的 影响 不 能 忽略 ,只 有 分 析 仿 生 鱼 在 静止 流 场 中 的 
运动 才能 单纯 地 研究 仿生 鱼 在 游 动 过 程 中 鱼 体 和 流 
体 之 间 的 作用 过 程 。 为 了 更 直观 地 了 解 仿生 鱼 在 流 
场 中 的 自主 游 动 行为 ,同时 也 为 了 验证 计算 结果 的 
合理 性 ,假定 鱼 体 处 于 一 静止 水 中 ,由 静止 开始 按照 
式 (9) 所 给 规律 在 水 中 运动 ,可 得 到 仿生 鱼 自主 游 
动 过 程 中 的 受 力 以 及 运动 情况 ,压力 场 分 布 、 涡 量 分 
布 以 及 等 值 面 图 如 图 3 ~ 图 5 所 示 。 


Pressure 


(a)t=3s 


(b)1=3.5s 


(0 1=4s 
图 3 不 同时 刻 压 力 分 布 
Fig.3 Pressure distribution at different moments 

从 图 3 可 以 看 出 ,仿生 鱼 在 自主 游 动 过 程 中 , 头 
部 位 置 受 到 的 阻力 较 大 ,在 所 有 时 刻 都 处 于 高 压 状 
态 , 鱼 体 两 侧 在 不 同时 刻 会 依据 鱼 体 摆动 状态 交替 
出 现 高 低压 区 域 。 仿 生 鱼 游 动 过 程 中 靠近 尾部 位 置 
在 摆动 方向 会 出 现 压 力 较 高 的 区 域 ,而 在 相反 方向 
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则 为 低压 区 域 , 这 些 压 力 差 的 存在 为 仿生 鱼 的 游 动 
提供 了 动力 。 同 时 在 尾 流 区 域 会 出 现 比较 明显 的 压 
力 扰动 区 域 和 洲 涡 脱离 现象 ,这 些 现象 与 相关 文献 
中 描述 鱼 的 运动 规律 比较 吻合 ,验证 了 本 人 研究 所 开发 
程序 的 正确 性 。 此 外 ,通过 仿生 鱼 自主 游 动 时 的 等 值 
面 分 布 图 可 以 看 出 ,仿生 鱼 在 自主 游 动 时 其 压力 和 涡 
量 分 布 是 三 维 的 ,压力 和 涡 的 分 布 在 尾部 两 侧 比 较 明 
显 , 而 在 仿生 鱼 中 心 线 附近 的 分 布 并 不 明显 。 因 此 ， 
采用 三 维 模拟 可 以 对 仿生 鱼 在 自主 游 动 过 程 中 的 受 
力 以 及 周围 流 场 的 变化 进行 较为 准确 的 描述 。 


Vorticity Magn 


(b)t=3.5s 
tude 


(c)t=4s 
图 4 不 同时 刻 涡 量 分 布 


Fig.4 Vorticity distribution at different moments 


2.2 来 流 方向 的 影响 


仿生 鱼 在 水 中 游 动 ,水 流 的 速度 和 方向 并 不 是 
保持 不 变 的 ,而 这 种 流 场 的 变化 必然 会 对 仿生 鱼 的 
运动 产生 一 定 的 影响 。 因 此 ,为 了 分 析 来 流 参 数 对 
仿生 鱼 自主 游 动 行为 的 影响 ,人 研究 了 在 来 流速 度 为 
0.1 m/s,0 分 别 为 90°、165°* 和 180°? 时 仿生 鱼 的 自主 
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游 动 和 受 力 情况 。 在 4s 时 各 9 角 下 的 压力 场 和 涡 
场 分 布 如 图 6 ~ 图 7 所 示 。 


Pressure 
25 
20 
IS 
10 
三 


-人 
© 

WO -10 
6 2 、 居 
-25 
化 3 
-40 
-45 


(a) 压力 分 布 
Vorticity Magnitude 


(b) 涡 量 分 布 
习 5 4=4s 时 等 值 面 图 


Fig.5 Jso-surface at 1 =4s 


Pressure 


(a) 0=90° 


(b) 0=165° 


(oj 0=180° 
图 6 流速 w=0.1 m/s 不 同 96 角 4s 时 压力 场 
Fig.6 Pressure distribution at 4 s with 
different 0 when wu =0.1 m/s 
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本 图 7 流速 w=0.1 m/s 不 同 9 角 4s 时 涡 量 分 布 

5 Fig.7 Vorticity distribution at 4 s with different 

CC 0 angels when u =0.1 m/s 

己 从 图 6 和 图 7 可 以 看 出 ,在 流体 速度 为 
0. 时 ns ,来 流 方向 垂直 于 鱼 自主 游 动 方向 时 ,在 仿 
生 鱼 的 周围 有 比较 明显 的 压力 变化 区 域 , 负 压 值 也 
较 大 。 当 来 流 方向 和 仿生 鱼 自主 游 动 方向 夹 角 变 为 
165° 时 , 鱼 体 周围 出 现 正 负 压 力 区 域 ,并 存在 一 定 的 
压力 梯度 。 当 来 流 方向 和 和 鱼 体 自主 游 动 方向 相反 
时 ,在 仿生 鱼 的 后 部 和 周围 可 以 出 现 比较 明显 的 压 
力 变化 区 域 ,而且 压 力 梯度 变化 范围 也 是 最 大 的 。 
另外 ,通过 计算 发 现 : 当 0=90" 时 在 4s 时 刻 仿生 鱼 
的 游 动 速度 为 -0.174 m/s; 当 09=165° 时 在 4s 时 刻 
仿生 鱼 游 动 速度 为 -0.1533 m/s;0=180°? 时 在 4s 
时 刻 仿 生 鱼 自主 游 动 速度 大 小 为 -0.187 m/s。 虽 
然 6=180° 时 流体 速度 分 量 全 部 作用 在 鱼 体 运动 方 
向 上 ,但 是 仿生 鱼 的 游 动 速度 却 是 最 大 的 。0 = 90° 
时 , 鱼 体 只 受到 垂直 方向 的 流体 分 量 的 影响 ,自主 游 
动 方向 的 作用 为 0, 但 是 自主 游 动 的 速度 却 小 于 9 = 
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180° 时 的 速度 。 而 当 0 =165° 时 ,仿生 鱼 在 游 动 方向 
和 垂直 方向 都 会 受到 流体 流动 分 量 的 作用 ,此 时 自 
主 游 动 速度 反而 要 比 0 = 90" 时 的 还 要 小 。 可 见 , 当 
流动 方向 和 鱼 体 运 动 方向 平行 时 ,可 以 获得 较 大 的 
自主 游 动 速度 ,运动 阻力 最 小 , 当 来 流 方向 与 鱼 体 游 
动 方程 垂直 (0 = 90°) 时 , 自主 游 动 速度 相对 0 = 
180° 时 要 小 ,而 游 动 方 向 与 来 流 方向 存在 一 定 夹 角 
(90° < 9 <180°) 时 ,自主 游 动 的 速度 最 小 ,阻力 
最 大 。 

在 同一 时 刻 不 同 9 角 时 , 鱼 体 不 同方 向 受 力 情 
况 如 图 8 ~ 图 9 所 示 。 
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图 8 不 同 6 角 4s 时 刻 鱼 体 表面 «方向 压力 分 布 


Fig.8 Pressure distribution in x direction at 4s 


with different 0 angels 
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图 9 不 同 0 角 4s 时刻 鱼 体 表 面 y 方向 压力 分 布 
Fig.9 Pressure distribution in y direction 
at4s with different 0 angels 

从 图 9 可 以 看 出 , 当 9 =90" 时 仿生 鱼 沿 长 度 方 
向 均 出 现 正 负 压力 区 域 , 且 正 负 压 差 较 大 ,分 布 存 在 
一 定 的 不 对 称 性 ,这 种 正 负 压力 区 域 的 产生 正 是 由 
于 涡 的 脱离 产生 的 ,也 是 鱼 向 前 运动 的 动力 。 此 外 ， 
在 y 方 向 鱼 体 两 侧 压 力也 存在 明显 的 不 平衡 ,两 侧 
正 负 压 差 较 大 , 迎 着 来 流 方向 的 那 一 边 正 负 压力 都 
大 于 背 流 区 域 。 当 0 = 165 "和 0 = 180。 时 鱼 头 部 位 
有 明显 的 高 压 区 域 , 尾 部 均 存 在 正 负 压 力 区 域 且 不 
对 称 , 鱼 体 在 长 度 方 向 上 交替 出 现 正 负 压 力 区 域 ,说 
明 仿 生 鱼 在 游 动 方向 受到 比较 大 的 驱动 力 , 两 种 角 
度 下 鱼 体 表面 的 压力 分 布 基本 相同 ,但 是 在 9 = 180° 
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时 仿生 鱼 尾 部 的 正 负 压力 差 较 大 ,从 yy 方向 鱼 体 表 
面 压力 分 布 图 可 以 看 出 , 鱼 体 上 下 表面 的 压力 分 布 
对 称 , 有 和 较 大 正 负 压 差 存在 ,9 =180° 时 压力 差 更 大 
一 些 。 因 此 ,通过 鱼 体 表 面 受 力 情况 可 以 看 出 , 鱼 尾 
部 的 正 负 压力 的 不 平衡 是 引起 鱼 运 动 的 主要 原因 ， 
而 且 这 种 压 差 越 大 ,仿生 鱼 所 受到 的 驱动 力也 会 
越 大 。 


3 结 论 


本 研究 针对 仿生 鱼 在 逆流 条 件 下 的 自主 游 动 行 
为 进行 仿真 模拟 ,人 研究 来 流 角度 对 仿生 鱼 游 动 行为 
的 影响 。 与 以 往 针对 仿生 鱼 游 动机 理 和 分 析 方 法 的 
研究 相 比 ,更 注重 考虑 外 界 流体 环境 对 仿生 鱼 游 动 
行为 的 影响 ,对 仿生 鱼 游 动 控制 策略 的 制定 具有 一 
定 的 参考 价值 。 通 过 计算 可 以 看 出 ,仿生 鱼 在 水 中 
自主 游 动 的 过 程 中 , 鱼 体 表面 会 交替 出 现 正 负 压 力 
区 域 ,其 尾 流 区 域 会 出 现 比 较 明显 的 涡 脱离 现象 ,这 
种 涡 的 脱离 会 给 鱼 提 供 前 进 的 动力 。 从 不 同 来 流 方 
向 时 的 仿生 鱼 的 运动 可 以 看 出 , 当 仿 生 鱼 游 动 方 向 
与 来 流 方向 平行 时 可 以 获得 比较 大 的 自主 游 动 速 
度 , 当 来 流 方 向 和 仿生 鱼 游 动 方向 有 一 定 角度 时 , 虽 
然 在 游 动 方 向 上 流动 分 量 较 小 但 是 鱼 体 自主 游 动 的 
速度 反而 要 小 。 男 外 ,从 仿生 鱼 体 表面 压力 分 布 图 
可 以 看 出 ,仿生 鱼 尾 部 明显 正 负 压力 不 平衡 的 存在 
是 产生 驱动 力 的 原因 ,在 仿生 鱼 运 动 存在 逆向 流体 
分 量 时 ,仿生 鱼 的 运动 控制 可 以 考虑 迎 着 水 流 运 动 
或 者 沿 迁 直 于 水 流 方向 运动 ,这 样 可 以 提高 运动 
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